Hysteresis behavior in current-driven stationary resonance 
induced by nonlinearity in the coupled sine-Gordon equation 
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Recently novel current-driven resonant states characterized by the 7r-phase kinks were proposed 
in the coupled sine-Gordon equation. In these states hysteresis behavior is observed with respect 
to the application process of current, and such behavior is due to nonlinearity in the sine term. 
Varying strength of the sine term, there exists a critical strength for the hysteresis behavior and the 
amplitude of the sine term coincides with the applied current at the critical strength. 
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Introduction. Recently the coupled sine-Gordon equa- 
tion has been intensively studied numerically and ana- 
lytically as a model of THz electromagnetic wave emis- 
sion from intrinsic Josephson junctions (IJJs). [l|, Q In 
current-driven emission without external magnetic field, 
novel resonant states characterized by the 7r-phase kinks 
were proposed for fairly large surface impedance Z, 
and the present author showed [3| that such states were 
stationary for any Z with large enough current. 

In these resonant states hysteresis behavior is observed. 
Strong emission which clearly breaks the Ohm's law only 
takes place in the current-increasing process, and the law 
almost holds in the current-decreasing process. Such dy- 
namical behavior is due to nonlinearity of the system, 
namely the sine term in the present equation. The main 
question of the present Letter is how such nonlinearity 
affects on hysteresis behavior in the current-driven reso- 
nance. May infinitesimal nonlinearity cause such behav- 
ior, or may there exist a critical value of it? 

In order to resolve this question, we introduce an artifi- 
cial parameter 7 in the equation and verify the strength of 
nonlinearlity continuously. As the parameter 7 decreases 
from the value in the original equation (7 = 1), width 
of hysteresis also decreases, and there seems to exist a 
nonvanishing critical value 7c with vanishing hysteresis. 

Model and formulation. When the capacitive coupling 
is not taken into account, IJJs are described by the cou- 
pled sine-Gordon equation, [6| 

dl,iJi = (1 - CA(2)) {df,i;i + I3dt>i:i + sinV-/ - J') , (1) 

with the layer index I and the operator A*^^-' defined in 
A^^^X/ = Xi^i — 2X1 -\- Xi^i. Quantities are scaled as 



x/Xc 



t' = Wpi, 



J' ^J/J,; t^p = c/ (yi^Ac) , (2) 



with the penetration depth along the c axis Ac, the 
plasma frequency in each layer Wp, and the critical cur- 
rent Jc. Using material parameters of Bi2Sr2CaCu208 
given in Ref. [7|, we result in a large inductive coupling 
C = 4.4 X 10^, and = 10 and /3 = 0.02 are taken. 

Neglecting temperature fluctuations and assuming ho- 
mogeneity along the y axis, we have the two-dimensional 



formula ([T]). Following our previous study, [g] width of 
the junction is chosen as — 86/xm, and the periodic 
boundary condition (PBC) along the c axis is consid- 
ered. Then, plasma velocity of the stationary state au- 
tomatically coincides with that of light in IJJs, which 
corresponds to the case with infinite number of junctions, 
though actual number of junctions N still affects on phys- 
ical properties even in the PBC, especially on response 
to the in-plane magnetic field. 'S] Here we concentrate 
on the simplest case iV = 4, which takes spatial inhomo- 
geneity of the superconducting phase into account. 

Since direct evaluation of electromagnetic wave emis- 
sion from edges of a thin sample to vacuum is quite 
complicated, we use a simplified version of the dy- 
namical boundary condition, where effects outside of 
the sample are only included in the relation between dy- 
namical parts of the rescaled electric and magnetic fields, 
E'l — ^ZB'i, with rescaled quantities E'l and related 
with ipi as dt'4'i = E'l and d^/tpi = (1 — ^A(^^)i3;, respec- 
tively. Here we take Z = 30, which gives strong enough 
emission close to the optimal value. [5| The sample along 
the X axis is divided into 80 numerical grids, and cal- 
culations are based on the RADAU5 ODE solver. ll| In 



order to investigate effect of nonlinearity in the sine term, 
Eq. ([ij is slightly modified as 

dl,iJi = (1-CA(2)) {df,ibi+pdt'^i+^sm^i - J') , (3) 

and the parameter 7 is controlled hereafter. 

Hysteresis in the original equation. First, the I-V 
curve of the original equation ([T|) is displayed in Figs.IUJa) 
and[IJb). As long as the fundamental and first harmonic 
modes are observed, large hysteresis is quite apparent. 
In the current- increasing process (Fig. [Ija)), the current 
rapidly increases as the voltage approaches the value cor- 
responding to the cavity resonance point given by the ac 
Josephson relation (broken lines) , 



2Li 



1.14 n [mV], 



(4) 



with the number of nodes n(= 1: the fundamental mode) . 
In the vicinity of emission peaks, the voltage exceeds the 
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FIG. 1: I-V curve (solid line) for the (a) current-increasing 
and (b) current-decreasing processes. Dotted line stands for 
the curve in another figure, and the solid-dotted lines denote 
the voltages corresponding to the cavity resonance points (H)). 



value given by Eq. which assumes perfect cavity res- 
onance. In experiments about 10% of frequency looks 
tunable in a single resonant branch by varying the cur- 
which is consistent with this result. 



rent, 



12 



13| which is consistent with this result. When 
a much larger value of Z is taken, the range of voltage 
becomes much smaller Q during varying the current in 
the region of strong emission. 

Emission intensity corresponding to the situations for 
Figs.[Ua) and[Hb) is plotted versus current in Figs.jUJa) 
andHJb), respectively. Strong emission observed in the 
current-increasing process is quite reduced in the current- 
decreasing process. Especially in the fundamental mode, 
the maximum intensity in the current-decreasing process 
is almost one order smaller than that in the current- 
increasing process. In such a case emission is expected 
to be invisible when experimental noise is overloaded in 
the reverse (current-decreasing) process, which may rep- 
resent "irreversible" emission in experiments. 

Hysteresis in the equations with varying nonlinearity. 
Next, the I-V curve of the modified equation ([3]) is in- 
vestigated. Here we concentrate on hysteresis behavior 
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FIG. 2: Emission intensity versus current for the (a) current- 
increasing and (b) current-decreasing processes. Dotted line 
stands for the curve in another figure. Circles and squares 
represent the emission peak in each node and in the current- 
increasing and current-decreasing processes, respectively. 



around the emission peaks in the fundamental {n = 1) 
mode, and only observe the n = 1 curves in the current- 
increasing process near the upper edges and the n = 2 
curves in the current-decreasing process near the lower 
edges. In Fig. |31Ia), parameter dependence of the upper 
and lower edges is visualized by various symbols for a 
wide range of parameters between 7 = 0.8 and 0.1. As 7 
decreases, the width of hysteresis, namely the difference 
of the voltages at the both edges in the current- varying 
processes, becomes smaller and smaller. Hysteresis is ob- 
served up to 7 = 0.2, while it is invisible at 7 = 0.1. 

Then, more precise measurement is made between 
7 — 0.2 and 0.1 as shown in Fig. [2Jb). Although the 
upper edges of the n = 1 curve decrease monotonically, 
the lower edges of the n — 2 curve exhibit non-monotonic 
behavior. The lower edges decrease up to 7 = 0.17, in- 
crease up to 7 = 0.15, and decrease again for smaller 7. 
Hysteresis behavior is observed up to 7 = 0.16, and from 
7 = 0.15 current at the upper and lower edges is the 
same, or the current-varying process becomes reversible. 
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FIG. 3: I-V curve for various values of 7 for (a) 7 — 0.8 
to 0.1 and (b) 7 — 0.2 to 0.1 (expanded figure around tiie 
irreversible-reversible boundary). Edges of the n = 1 curve 
in the current-increasing process and the n = 2 curve in the 
current-decreasing process are visualized by various symbols. 



similar behavior may also be observed. Numerical study 
along this direction is now in progress. [l5| 

Frequency of resonance / in Eq. ^ is derived from the 
linearized version of Eq. ([T]), [3] and the role of nonlin- 
earity is to generate the 7r-phase kinks Q and hysteresis 
behavior. The present result shows that the hysteresis 
behavior is also nontrivial. We consider this behavior is 
general and can be regarded as a prototype of hysteresis 
in harmonic oscillations induced by nonlinearity. 

Summary. In the present Letter we numerically inves- 
tigate characteristic behavior of novel current-driven res- 
onant states in the coupled sine-Gordon equation. These 
states were proposed theoretically as strong emission 
states in THz electromagnetic wave emission from intrin- 
sic Josephson junctions, and such emission depends on 
the application process of current. That is, hysteresis is 
observed in emission behavior driven by applied current. 
Such irreversible behavior is due to nonlinearity of the 
equation, and we vary the strength of the sine term. 

As long as the fundamental resonant mode is observed, 
the maximum intensity of emission (and consequently 
the maximum value of the current) decreases monotoni- 
cally as the strength of nonlinearity decreases, and at the 
nonvanishing critical value hysteresis behavior disappears 
and the amplitude of the sine term coincides with that of 
the applied rescaled current at the emission peak. From a 
physical point of view, superconductivity and hysteresis 
behavior vanish at the same time. 
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Discussions. These results suggest that the boundary 
between the irreversible and reversible behavior locates 
between 7 = 0.16 and 0.15, but further precise calcula- 
tion in search for the next digit of 7 may not be produc- 
tive. Instead, we point out that the current at 7 = 0.15 is 
J' = 0.150 on the edges, or that the amplitude of the sine 
term coincides with the constant term at 7 = 0.15 in the 
equation ([3]). This fact strongly suggests that 7 — 0.15 
is the irreversible-reversible boundary. Quite recently re- 
versible THz wave emission from IJJs was reported, , 
which may be related with the present finding. 

Introduction of the parameter 7 can be regarded as 
modification of the critical current, namely Jc — jJc- 
Then, j — J' means that the current J is equal to the 
modified critical current jJc, where superconductivity 
breaks down and the present model is not justified any- 
more. It is natural that hysteresis behavior of the model 
may also change there. In the original equation ([1]), the 
situation J > Jc occurs in higher harmonic modes and 
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